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We study the massless scalar wave propagation in the time-dependent Schwarzschild black hole 
background. We find that the Kruskal coordinate is an appropriate framework to investigate the 
time- dependent spacetimes. A time-dependent scattering potential is derived by considering dynam- 
ical black hole with parameters changing with time. It is shown that in the quasi-normal ringing 
both the decay time-scale and oscillation are modified in the time-dependent background. 

I. INTRODUCTION 

The study of the Quasi-Normal Modes (QNM) developing outside black holes plays an important role in black hole 
physics and astrophysics. In virtue of previous works, we now have the schematic picture regarding the dynamics of 
waves outside a black hole. A static observer outside the black hole can indicate three successive stages of the wave 
evolution. The first stage is the exact shape of the wave front, which depends on the initial pulse. This stage is 
followed by the quasi-normal ringing, which has characteristic frequencies directly connected to the parameters of the 
black hole. This quasi-normal ringing is regarded as a direct evidence of the existence of black holes and is expected 
to be detectable by future gravitational wave detectors. Finally, at late time, the quasi-normal ringing is replaced by 
a relaxation process, which is known as the wave tails. Wave tail is of interest in its own right since it has various 
applications from the no-hair theorem to the mass inflation scenario. 

Quasi-normal modes of massless scalar, electromagnetic, and gravitational fields in asymptotically flat spacetimes 
have been studied for over thirty years. For a review, see [1]. It was shown that the quasi-normal ringing decays 
exponentially with time. The oscillation period and the decay time scale are associated with the real part and 
the imaginary part respectively of the so called complex quasi- normal frequency [2]. The tail phenomenon was first 
investigated by Price [3] for both Schwarzschild and Reissner-Nordstrom black hole. It was shown that for an observer 
at a fixed position the field dies off with a power-law tail. This behavior was confirmed by numerical studies [4] . The 
evolution of massive scalar field in the asymptotically flat black hole background is also of interest. Recently this 
topic has attracted a lot of attention [5]. 

The study of wave dynamics is not limited in the asymptotically flat spacetimes. Motivated by the inflation, the 
quasi-normal modes in de Sitter (dS) space were studied [6, 7]. One of the most significant differences between the 
QNMs in the dS space and the QNMs in the asymptotically flat space is that in the dS space the oscillation always 
decays exponentially, the power-law tail is not in presence. Inspired by the discovery of the AdS/CFT correspondence, 
the investigation of the quasi-normal modes in the Anti de Sitter (AdS) black hole has been appeahng [8, 9, 10, 11] 
recently. It was shown that the decay of the perturbation is always exponential for the AdS black holes. 

Another intriguing aspect of studying QNM has been suggested recently. It was argued that QNM could be a probe 
to investigate the quantum properties of black holes. The first speculation was put forward by Hod who pointed out 
that the QNMs of a Schwarzschild black hole conform to the existence of the Bekenstein-Hawking temperature and 
imply emission in terms of a fixed area quantum [12]. This study was subsequently extended to other black holes such 
as the d-dimensional spherically symmetric black holes [13], near extreme Schwarzschild-de Sitter and Kerr black holes 
[14], etc [15]. Though this kind of consideration is still somewhat speculative, it possibly relates classical, macroscopic 
properties to quantum information of black holes. 

All these previous works on quasi-normal modes have so far been restricted to time-independent black hole back- 
grounds. It should be realized that, for a realistic model, the black hole parameters change with time. A Schwarzschild 
black hole gaining or losing mass via absorption or evaporation is a good example. The more intriguing investigation 
of the black hole QNM calls for a systematic analysis of time-dependent spacetimes. Recently the late time tails under 
the influence of a time-dependent scattering potential has been explored in [16], where the tail structure was found 
modified due to the temporal dependence of the potential. The motivation of the present paper is to explore the 
modification to the QNM in the time-dependent spacetimes. Instead of putting an effective time-dependent scattering 
potential by hand as done in [16], we will introduce the time-dependent potential in a natural way by considering 
dynamic black holes with black hole parameters changing with time due to absorption and evaporation processes. We 



2 



will study numerically the temporal evolution of massless scalar field perturbation, especially the QNM in different 
time-dependent situations. 

The outline of this paper is as follows. In Sec. II we first go over the conventional treatment for numerical study of 
the wave propagation in the stationary black hole background. We will point out the difficulty if it is to be applied to 
the time-dependent cases. Then we will introduce a new treatment by employing the Kruskal coordinates which we 
find to be advantageous. In Sec. Ill, we present the numerical study of the QNM in the time-dependent background. 
Sec. IV is a brief summary and discussion. 

II. FIELD EVOLUTION IN TIME-INDEPENDENT BACKGROUND 
A. The conventional treatment 

In this section, we first present a brief review on the conventional treatment of the evolution of massless scalar field 

propagating in a stationary Schwarzschild black hole background. For more details, please refer to [4]. 

The metric of a stationary (time-independent) Schwarzschild black hole in the general coordinates is given by 

ds^ = -fdf + f-^dr^ + r'^dn'^, (1) 

where / = 1 — 2M/ r. The propagation of waves in curved spacetimes is governed by the Klein-Gordon equation 

= 0. (2) 

Substituting the metric, Eq. (2) becomes 

^$,ft - (rV,r + 2r/) - rV^.rr - cot e^,g - ^,00 - -4^^,00 = 0. (3) 
/ sm 

Since the background is spherically symmetric, each multipole of a perturbation field evolves separately. Hence one 
can write 

$(i,r,0,</>) =^((f,r)y,„(0,</>)/r, (4) 
and obtains a wave equation for each multipole moment 

'1(1 + 1) , /,r 



/ V,M - f->P,rr - f,r->P,r 

Using the tortoise coordinate r* defined by 



the wave equation becomes 
where the effective potential 

Eq. (7) can be recast as 
where the null coordinates 



■ip = 0. (5) 

dr* = f-^dr, (6) 
ttt - tr'r^ + = 0> (7) 



Vir) = f 



(8) 



i^,uv + \v{r)i; = 0, (9) 



u = t-r*,v = t + r*. (10) 
To write Eq. (9) into the discrete form, one can construct a functional 

IbP]= jjdudvF[u,v,il),il)^u,i>,v], (11) 
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and the differential equation Eq. (9) is recovered from the functional variation 



Sip 



0, (12) 



or 



if the function F takes the form 



d_fdF_\ ^f_^\^^_r. 



F[u,v,il},il},u,'ipA = i',utp,v - ^V{r)ip'^. (14) 

Using the three point formula for the first-order derivatives and the trapezoid rule, one can convert Eq. (11) into 
a summation form 



I[4)] = AuAv ^ 



AuAv 



(15) 



where ijjij denotes ip{ui,Vj). 

Now each discrete ipij is treated as an independent variable, and the functional variation 57[^]/5^ = becomes a 
partial derivative 



Substituting Eq. (15), Eq. (16) becomes 



= 0. (16) 



ipi-i,j+i + A+ij-i - A-ij-i - A+ij+i - ^AuAvV{r)'4)ij = 0. (17) 



Taking the following notations 



i>E = 1pi+l,j-l, tpW = V'i-lj + l, 



(18) 



and replacing the midway quantity between two neighboring points by the average of two lattice points 



one arrives at 



V;;v = V'M- + V'E - V's - AuAvV{r) + = o. (20) 

8 

It is straightforward to setup a rectangular grid of u and v. Starting from initial data on u = uq and v = vq, 
integration can proceed to the northeast, as described in [4] and numerical result of Eq.(9) can be obtained. 

The conventional treatment is powerful in studying the wave propagation in stationary black hole backgrounds, 
however it is difficult to be extended to investigate the time-dependent spacetimes. For a dynamic Schwarzschild 
black hole, the mass of the black hole is a function of time M = M{t), then we get another term —{f,t/.P)iJ-',i in Eq. 
(5) for each multipole moment. Due to the appearance of the first-order derivative associated with time, the neat 
expression Eq. (7) and Eq. (9) are impossible to be obtained. Such an equation is difficult to handle. Besides the 
tortoise coordinate r* now is not only a function of r, but also a function of t. These difficulties make the conventional 
numerical treatment no longer appropriate in the time dependent background. 
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B. Wave equation in the Kruskal coordinates 



To avoid these difficulties, we use the Kruskal coordinates instead of the general coordinates Eq. (1). As we 
know, the metric terms of the timc-likc coordinate r and the space-like coordinate p arc the same in the Kruskal 
coordinates, hence the mass of the black hole M{t) only appears in the angular part of the wave equation. Thus in 
the time-dependent case, the investigation will be easier. This will be exhibited in Sec. III. Here we first present the 
investigation in the stationary black hole background in order to test the validity of applying the Kruskal coordinates 
in studying the wave propagation. 

In the Kruskal coordinates, the metric for the Schwarzschild black hole is given by 

= 6-/2^(^^2 _ ^^^2) _ ,.2^f^2_ (21) 

r 

where r and p are the time-like and space-like coordinates respectively. In the region r > 2M, they are defined by 



r _ ^ ^,/4M gij^i^ ±_ p= lL_ I ^r/AM ^osh -L. (22) 

2M 4M y 2M 4M ^ ' 

The propagation of the massless scalar field is governed by the wave equation 

2r sin 9r^r^,T + t'^ sin ^$,rr - 2r sin Or^p^^p - sin B^,pp 

_ cos Qe-^l^^^ e - ^ sin^^e-'■/2^$ - ^g-''/^^* = Q. (23) 

r ' r r sm & 

Since the background is spherically symmetric, we can write ^{t, p,0,(j)) = ^^{t, p)Yim{0,(l))/r, and we get the wave 
equation for each multipole moment 



32M^l{l + l)e-''/2^ 



^3 1^ lyt 



V' = 0. (24) 



Analogous to the null coordinates used in [4], we make the following variable transformations 

u = T — p, V = T + p, (25) 



and the wave equation becomes 



where 



i^,uv + V^ = 0, (26) 



The expression of r,„„ is simple when M is a constant. For the stationary black hole, we have 
It is straightforward to write Eq. (26) into the discrete form 



(28) 



i'N = i'w+^E-i's- AuAvV^I^^^^^^. (29) 

The only difficulty of developing the code is to retrieve r and t from Eq. (22) and Eq. (25) for each point on 
the grid. In a typical setup, the approximate range of u and v arc from —10*^ to —10"^" and from 10^ to 10^° 
respectively. In order to finish the integration over such a big region in a reasonable duration of time while sticking 
to an acceptable precision, we divide the whole grid into small blocks. The step length Au and Av vary from block to 
block in accordance with the requirement. Starting from the southwest, the small blocks are resolved one by one to 
the northeast, the output of each resolved block becomes the input of the forthcoming blocks. We use 'ip{v = vq) = 
and a Gaussian pulse on u = uq as the initial data for all of our calculations. 

The numerical evolution of the / = 2 mode is displayed in Fig. 1. For comparison, we also exhibit the curve 
obtained in the general coordinates using the conventional treatment. Since the amplitude of the oscillations is 
physically irrelevant, we concern ourselves on the quasi-normal frequencies of the oscillations. The real part of the 
quasi-normal frequency ojr and the imaginary part of the quasi- normal frequency oji read from Fig.l and the quasi- 
normal frequencies of a few other modes are listed in Table I. As one can see from the table, numerical result got by 
employing the Kruskal coordinates agrees well to that of the conventional treatment. This shows that our treatment 
using the Kruskal coordinate is valid to investigate the wave propagation in the black hole background. 
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FIG. 1: Temporal evolution of the field tp ai a. fixed radius r = 2. The mass of the black hole M = 0.5 and the multipole 
index / = 2. The upper curve is obtained from calculation done in the Kruskal coordinates, the lower curve is calculated in the 
general coordinates. 



TABLE I: Quasi-normal frequencies for the massless scalar perturbation in the stationary Schwarzschild background. The mass 
of the black hole M = 0.5. 



General Coordinates 

(JJR OJl 

0.586" -0.1956'' 

0. 968 -0.1932 
1.351 -0.1926 

1. r:i() -0.11)2 i 



Kruskal Coordinates 

0.587 -0.1947 

0.968 -0.1932 

1.355 -0.1926 

1.735 -0.1922 



m is 0.56 in [4]. 
'It is -0.19 in [4]. 



III. FIELD EVOLUTION IN TIME-DEPENDENT BACKGROUNDS 



If the mass of the black hole has a general dependence on t, then M is a function of r and p in the Kruskal 
coordinate. Since the metric gTj- = Qpp, the mass M only appears in the angular part of the wave equation. This is a 
big advantage of using the Kruskal coordinates in study of the time-dependent case. The wave equation has the same 
form as Eq. (26), however now the second-order derivative r,„„ is replaced by 

- ^ + 7^ • (30) 

The derivatives M,„, M,^, and M^uv are given by 

M 

M,u = ^{t,r-t,p), (31a) 
M 

Mv = Yii,r+t,p), (31b) 

M , y y , M , , , , 

M^uv = -r{i,r - i,p) + ^(i.rr - i,pp), (31c) 



6 




10 20 30 40 



80 90 100 



FIG. 2: Temporal evolution of the field in the background of a Schwarzschild black hole for I — 2 and r = 2. The mass of the 
black hole M{t) = Mo ± at, where Mo = 0.5 and a = 5 x 10~* are constant coefficients. The field evolution for M = Mo + at 
and M = Mo — at are shown as the top curve and the bottom curve respectively. For comparison, the oscillations for M = Mo 
are also displayed (the middle curve). 



where M and M are the first-order and second-order derivatives with respect to t. The terms t^r, t,p, t,TT^ and t^pp 
can be obtained directly from Eq. (22) 



^ _ 4M2cosh^^ 
p{M-Mt) ' 



T{M-Mt) 

UxiAi^iM - Ml)- + 2M-Mt + iMM(M - Mt) 



2M2(M - Mt) 
coth 4|y(M - Mtf + 2M^Mt + AMM{M - Mt) 



2M'^{M-Mt) 



(32a) 
(32b) 

(32c) 

(32d) 



We now present the result of our numerical calculations in the time-dependent black hole background. In the first 
scries of simulations, wc consider the simple situation by choosing the mass of the black hole M = Mq ± at, where 
Mq and a are constant coefficients. The results are shown in Fig. 2, 3, and 4. The modification to the QNM due 
to the time-dependent background is clear. When M increases linearly with respect to t, the decay becomes slower 
compared to the stationary case, which corresponds to say that the uji is not a constant and decreases with respect 
to t in the time-dependent situation. The oscillation period is no longer a constant as that for the stationary black 
hole, it becomes longer with the increase of time. In other words, the real part of the quasi-normal frequency ur 
decreases with the increase of time. When M decreases linearly with respect to t, compared to the situation in the 
stationary black hole, we have observed that the decay becomes faster and the oscillation period becomes shorter 
which corresponds to both loj and cor increase with the increase of time. 

We have also extended our discussion to a more realistic model, an evaporating Schwarzschild black hole with the 
mass determined by [17] 



M 



where ao is a constant coefficient that takes different values for black holes with mass M 3> 10"'^^, 
mass M <C lO^^g. Since the wave equation Eq. (23) is invariant if one makes the rescale M 



(33) 

and those with 
aM, r ar, and 
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FIG. 3: The same plot as Fig. 2. For clarity, the connected maxima of the oscillations are displayed. The imaginary frequency 
uji can be read from the decay time scale which is proportional to l/oij. 
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FIG. 4: The frequency u>r determined from the oscillation period in Fig. 2 is shown as a function of t. 



t — > at, the value of ag is not very important here. One can obtain the mass of the black hole as a function of t from 
Eq. (33) as 



M{t) = [3ao(6-0]'^'. 



(34) 



where h is an arbitrary constant. Results of the numerical calculations are shown in Fig. 5, 6, and 7. Different 
from the stationary black hole case, for the evaporating black hole both real and imaginary parts of quasi-normal 
frequencies and ojr increase with respect to t, in consistent with the simple case above. 



8 



-10 - 



-15 - 



stationary black hole 
evaporating black hole 




10 



20 



30 



40 



50 



f 



FIG. 5: Temporal evolution of the field in the background of a evaporating Schwarzschild black hole for I = 2 and r = 2 
(lower curve) . We use ao = 2.011 X IQ-"* [17] and b = 207.2 in our calculation. For comparison, the I = 2 mode in the stationary 
background with M = 0.5 is also displayed (upper curve). 
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FIG. 6: Same as Fig. 5. Only the connected maxima of the oscillations are displayed. The imaginary frequency uJi can be read 
from the decay time scale which is proportional to l/cJi. 



IV. SUMMARY AND DISCUSSION 



Wc have studied the evolution of the massless scalar field in the time-dependent Schwarzschild black hole back- 
ground. We have found that the Kruskal coordinates is an appropriate framework to investigate the wave propagation 
in the time-dependent spacetimes. Despite imposing the effective time-dependent scattering potential in the wave 
equation by hand, in our study we have tried to derive the time-dependent potential in a natural way by considering 
dynamic black holes with black hole parameters changing with time. In our numerical study, we have found the 
modification to the QNM due to the temporal dependence of the black hole spacetimes. The decay and oscillation 
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FIG. 7: The frequency lur determined from the oscillation period in Fig. 5 is shown as a function of t for the evaporating black 
hole. For the stationary black hole ujr is a constant. 
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FIG. 8: The change of the effective potential with time for stationary black hole and dynamical black holes. 

timescale are no longer constants with the evolution of time as that in the stationary black hole case. In the absorption 
process, when the black hole mass becomes bigger, both the real and imaginary parts of the quasi-normal frequencies 
decrease with the increase of t. However, in the evaporating process, when the black hole loses mass, both the real 
and imaginary parts of the quasi-normal frequencies increase with the increase of t. 

For the black hole evaporating case, black hole mass decreases with the increase of time and so does the effective 
potential, see Fig. 8. The correspondence between the decrease of decay time scale and decrease of the effective 
potential with the evolution of time is qualitatively in consistent with the tail behavior exhibited in [16], where 
V{x, t) = l/{x°'t^) and the power-law indices found to be —4.04 and —5.08 for /? = 1 and (3 — 2, respectively, showing 
that smaller potential (with bigger /3 for the same t and fixed x) leading to faster decay (bigger absolute value of 
power-law index). 



evaporating black hole 
stationary black hole 
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